Atom counting in ultra-cold gases using photoionisation and ion detection 
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We analyse photoionisation and ion detection as a means of accurately counting ultra-cold atoms. 
We show that it is possible to count clouds containing many thousands of atoms with accuracies 
better than 7V _1//2 with current technology. This allows the direct probing of sub-Poissonian number 
statistics of atomic samples. The scheme can also be used for efficient single atom detection with 
high spatio-temporal resolution. All aspects of a realistic detection scheme are considered, and we 
discuss experimental situations in which such a scheme could be implemented. 

PACS numbers: 03.75.-b, 03.75.Be, 32.80.Fb, 39.90.+d 
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I. INTRODUCTION 

The term "atom optics" arises from the analogy be- 
tween experiments exhibiting the wave-like properties of 
light with those demonstrating the wave-like properties 
of matter |l|. Due to the vanishingly small de Broglie 
wavelength of room temperature atoms it is necessary 
to cool gases to microKelvin temperatures before these 
properties are observable. Before 1995, experiments in 
atom optics used clouds of laser-cooled thermal atoms 
that can be considered similar to thermal light sources in 
optics, with coherence lengths of less than a micron. 

The realisation of Bose-Einstein condensation in 1995 
HI 01 El provided the first intense sources of coherent mat- 
ter waves, analogous to the invention of the optical laser. 
Many of the early experiments with BECs concentrated 
on their wave-like properties, which in general can be well 
described by the Gross-Pitaevskii equation for the clas- 
sical field ip( x ) s» (^>(x)). Some notable experiments in- 
clude the interference of independently prepared conden- 
sates and the demonstration of atom lasers 0, S H 13 • 
These experiments rely on the first-order coherence of the 
condensate, and in some sense can be considered to be 
coherent "classical atom optics" . 

Recently there have been an increasing number of ex- 
periments focusing on the generation and measurement of 
non-trival higher order coherences in matter waves that 
are beyond the scope of Gross-Pitaevskii mean field the- 
ory. These experiments could be said to be in the area 
of "quantum atom optics" E3 • The development of this 
new area bears a close resemblance to the development 
of quantum optics following the invention of the optical 
laser. Just as quantum optics relies on single photon de- 
tectors, quantum atom optics will utilise detectors with 
single atom resolution, such as the technique we describe 
in this paper. 

One of the earliest experiments probing higher-order 
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coherence in BECs was in the measurement of three-body 
loss from a BEC, from which the correlation function 
</ 3 )(x) could be inferred E2- There have been several 
measurements made of the local second-order correlation 
function g( 2 '(x.) for a number of systems, including a ID 
strongly interacting Bose gas I12lll3l , a quasi-condensate 
[l4j . and a thermal Bose gas |15|. The value of g^ 2 \x) 
for an atom laser has been measured in ref. |l6j| . 

There has also been strong interest in relative number 
squeezing. This is where a state is generated in which the 
number difference between two or more atomic samples is 
sub-Poissonian. Javanainen and Ivanov 

El 

showed that 

splitting condensates by adiabatically turning up a tun- 
nel barrier can produce number squeezed atomic states. 
Experimentally, Orzel et al. loaded a BEC into a 
one-dimensional optical lattice and observed the degra- 
dation of interference fringes. This was interpreted to be 
caused by increased phase-fluctuations due to reduced 
number fluctuations at the lattice sites. A similar ob- 
servation was made in the first demonstration of a 3D 
Mott insulator state Il3l and further elucidated by other 
related experiments [2Cj . 

Correlations and relative number squeezing have been 
predicted in collisions and four-wave mixing in conden- 
sates 01 ; as wen as m downconversion of a 
molecular BEC in both the spontaneous [25| and stim- 
ulated [2^| regimes. While experiments have observed 
four- wave mixing processes (2^, and matter- wave am- 
plification p9|. there has been no direct measurement of 
these correlations. 

Experimentally, other high-order correlations have 
been directly detected by suitable analysis of absorp- 
tion images from an ensemble of experiments [3tt l3l| . 
Other schemes such as the "quantum tweezer" [33 can 
deterministically extract sub-Poissonian atomic samples 
from a BEC. Optical dipole traps along with fluorescence 
imaging have been used to demonstrate sub-Poissonian 
statistics in small samples of atoms |33| . 

In order to directly probe atom number statistics of 
cold quantum gases, one would ideally like an atom 
counter with accuracy at the single atom level. Exist- 
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ing single atom detection schemes can be divided into 
two categories: optical and contact. One of the sim- 
plest optical schemes is based on resonance fluorescence 
detection [33L 133. l35| . Assuming negligible background, 
the shot noise limited atom number uncertainty in a flu- 
orescence measurement scales as y/N/ayr, where N is 
the number of atoms, a is the collection efficiency of the 
detector, 7 is the mean scattering rate and r is the inte- 
gration time |36| . With a long enough integration time it 
is possible to measure atom numbers with sub-Poissonian 
precision. However, limitations such as atom heating re- 
strict the capability of sub-Poissonian fluorescence detec- 
tion to atom numbers of order < fOO |33|. Absorption 
measurements offers more favourable scaling to higher 
atom numbers |36| but achieving good absolute accuracy 
remains technically difficult. High finesse optical cavi- 
ties offer excellent sensitivity for single atom detection 
[37L l38l |3t| and can be used for slow counting of a large 
number of atoms. Cavities offer an improvement in signal 
to noise of y/~T, where T is the cavity finesse, compared 
to single pass measurements for a fixed amount of heat- 
ing [I(J, |4y . Moderate cavities have also been shown to 
aid fluorescence and absorption measurements |3^ |. It 
is difficult however to apply cavity techniques to count- 
ing large atom numbers and long integration times are 
required. 

The second category of single atom detectors is based 
on contact methods. These generally involve a charged 
particle or metastable atom with high internal energy 
impacting on a charged particle detector. This initiates 
a cascade of electron emissions, which are amplified to 
produce a macroscopic current pulse. This method is 
particularly useful for detecting metastable He atoms |42| 
where the arrival statistics of atoms released from a BEC 
allowed for matter wave Hanbury-Brown Twiss correla- 
tions to be observed While this detection method 
is simple, metastable atoms remain difficult to condense 
and most cold atom experiments use alkali atoms. 

In this paper we describe an atom detection scheme 
based on photoionisation and ion detection. This scheme 
overcomes some of the difficulties associated with optical 
detection, can be applied to ground state alkali atoms and 
is capable of extending sub-Poissonian counting to much 
larger numbers of atoms. Photoionisation 0] and ioni- 
sation 44] have previously been considered as a means of 
efficient atom detection. However, these proposals have 
not included many factors present in realistic detectors. 
Firstly, we note that the best absolute detection efficien- 
cies of room temperature charged particle detectors are 
of order 0.8-0.9. (Superconducting tunnel junctions can 
in principle provide unity detection efficiency |45| but are 
impractical for many experimental setups.) This has se- 
rious consequences for any experiment which endeavours 
to count a large number of atoms, but as we show in this 
paper, provided the detection efficiency is greater than 
0.5, it is still possible to count atoms with accuracies 
better than \j\f~N . The achievable accuracy depends on 
several factors including the detector calibration, back- 
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FIG. 1: (a) In the proposed atom counting scheme atoms are 
ionised from an optical dipole trap via a three-photon process. 
The resulting ions are accelerated into an ion detector, (b) 
The ionisation scheme showing the relevant energy levels of 
87 Rb (not to scale). The scheme consists of 2 x 778 nm 
photon excitation to the 5D 5 / 2 state, followed by ionisation 
with a Nd:YAG laser, which also forms the dipole trap. 



ground count rates, detector pulse resolution times, ion- 
isation rates and the uncertainties in these parameters, 
all of which are included in our analysis. 

This paper is organised as follows: in Section II we de- 
scribe the experimental arrangement necessary for atom 
counting and how it can be applied to a cloud of trapped 
87 Rb atoms. Section III addresses the statistics associ- 
ated with the detection scheme along with the require- 
ments on the detector efficiency. In Section IV we dis- 
cuss detector calibration and describe numerical simula- 
tions of the calibration including the effects mentioned 
above. Section V looks at the application of the scheme 
in the presence of losses and we follow with a discussion 
of potential applications. Appendix lAlcontains a table of 
symbols and Appendices [B] and El detail the derivations 
of certain expressions used in the paper. 



II. SCHEME 

The method we propose consists of slowly photoion- 
ising a sample of atoms and accelerating the ions to a 
charged particle detector. Figure 1(a) shows the elements 
of the scheme. A cloud of cold atoms is held in an opti- 
cal dipole trap where it is illuminated by the ionisation 
lasers. Our discussion focusses specifically on the case of 
87 Rb atoms but the scheme could be modified to detect 
other alkali atoms by the correct choice of lasers. Figure 
1(b) shows the relevant energy levels and transitions for 
87 Rb. 

The ionisation scheme consists of two stages. The first 
is two-photon excitation with a low intensity 778 nm laser 
from the 5S!/ 2 state to the 5D 5 / 2 state. The second is 
ionisation from this state with a high intensity infra-red 
laser such as a Nd:YAG laser. The 5D 5 / 2 state has a 
photoionisation cross-section approximately three orders 
of magnitude higher than the ground state |4rj| so ioni- 
sation can readily be achieved with a cw laser of wave- 
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length shorter than ~ 1200 nm rather than with pulsed 
lasers 01 ■ I° n detection is achieved after appropriate 
ion collection optics using a standard ion detector such 
as a channel electron multiplier (CEM) or discrete dyn- 
ode detector. 

The atom counting scheme works best if the pulses pro- 
duced by ions striking the ion detector are well resolved 
in time. This means the ionisation rate iVj should be 
much lower than the inverse of the pulse resolution time 
of the detection system l/r r . In Section V we will see 
that this means that for typical detector efficiencies and 
trap lifetimes, and atom numbers of the order of 10 3 or 
10 4 , the ionisation should take place on a time scale of 
order several ms. In order to confine the atoms during 
this time, the Nd:YAG laser is focused to create a dipole 
trap. This localises the atoms in the highest intensity 
region of the beam. The intensity must be high enough 
to produce a trap sufficiently deep that no atoms are lost 
due to heating as the ionisation takes place. The ma- 
jor source of heating is the spontaneous decay of atoms 
excited to the 5D 5 / 2 state. A high Nd:YAG laser inten- 
sity thus helps eliminate trap loss not only by creating a 
deep trap, but also by ensuring that an atom excited to 
the 5D 5 / 2 state has a substantially higher probability of 
being ionised than of decaying back to the ground state. 
Atoms that decay from the 5D 5 / 2 state can end up in 
several ground state sublevels. A benefit of the dipole 
trap is that it traps all of these sublevels, in contrast to 
magnetic fields, which can only trap a restricted class of 
sublevels. 

The rate of ionisation can be controlled with the 778 
nm laser intensity. In Section V we show that for a given 
atom number and trap loss rate there is an optimal ioni- 
sation rate that minimises the uncertainty in the inferred 
atom number. In order to hold the ionisation rate ap- 
proximately constant at the ideal value, the 778 nm laser 
intensity can be exponentially ramped to compensate for 
the reduction in trapped atom number as ions are pro- 
duced. However we also show that the scheme works well 
with constant laser intensities. 

This state selective scheme also represents an efficient 
single atom detector The ionisation lasers can be 

tightly focussed to detect atoms at a well defined spa- 
tial location. The two photon excitation rate is propor- 
tional to the squared intensity of the 778 nm laser so 
that scattered light is less likely to heat atoms not in the 
detection region. Additionally, for single atom detection 
the spatial intensity profiles of the two lasers can be cho- 
sen so that the repulsive optical dipole potential of the 
blue-detuned 778 nm laser can be compensated for by 
the attractive potential of the high power red-detuned 
ionising laser. This leads to minimal perturbation of un- 
detected atoms which may for instance be magnetically 
detuned or trapped in a different internal state. These 
factors, and the relatively quick ionisation times, make 
this a flexible scheme well suited to single atom detection 
with excellent spatio-temporal resolution. 



III. COUNTING WITH IMPERFECT 
DETECTORS 

We now consider how accurately an atom number can 
be determined using a scheme in which ions are counted 
with a non-unity efficiency detector. The overall ion de- 
tection efficiency is 

Vi = VdetPr, (1) 

where rjdet is the efficiency of the ion detection system 
and is the product of the ion detector efficiency rj' det and 
the collection efficiency rj co u . The quantity p r is the prob- 
ability that an ion is temporally resolvable from the other 
ions detected, and is ionisation rate dependent. 

If JVj ions are produced, the number of ions counted, 
Ni, is described by a binomial distribution. Provided 
r}i(l — rji)Nj is greater than ~ 5, this is well approximated 
by a normal distribution with mean rjiNj and standard 
deviation y/r}i(l — r]i)Ni |49j. In a given experiment, the 
number of atoms ionised can be inferred to be 

N Lmf = —, (2) 
Vi 

with an uncertainty 

\ N t {l- ni ) Nfal 

hinl ~ V ~oT' () 

where a m is the uncertainty in rji. 

In order to compare this uncertainty to Poissonian fluc- 
tuations, we define the normalised count uncertainty: 



The quantity K is related to the measurement Fano 
factor discussed in Section IVll but includes an additional 
systematic contribution due to the uncertainty in the de- 
tection efficiency (the second term under the square root 
in Eq. |3J). k = 1 implies the count uncertainty is equal 
to the uncertainty inherent in Poissonian statistics. In 
order to count ions with error less than y/Ni, n must be 
less than one. 

In Fig. |21 contours of the normalised count uncertainty 
K are plotted as a function of r\i and <r Vi for Ni = 10 3 . 
From this figure it can be seen that for a detection sys- 
tem in which the detection efficiency is known exactly, 
i.e. a rji = 0, the overall ion detection efficiency rji need 
only be greater that 50% to achieve k less than 1. As 
the uncertainty in the detection efficiency increases, the 
required detection efficiency to achieve K < 1 also in- 
creases. Thus, to achieve a low value of n, we require 
a high efficiency detector that is well calibrated. In the 
next section we examine how this can be done. 

IV. DETECTOR CALIBRATION 

To infer the number of atoms ionised and its uncer- 
tainty from an ion count, we need to know the overall ion 
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FIG. 2: Contours of the normalised count uncertainty k for an 
atom counting system with ion detection efficiency rji ± a Vi , 
for ion number Ni — 10 3 . 



detection efficiency rji and its uncertainty <r ni . From Eq. 
© this means we need to know the ion detection system 
efficiency rjdet and the ion temporal resolution probabil- 
ity p r , each with their respective uncertainties. We first 
consider how to measure rjdet to high accuracy, or in other 
words, how to calibrate the ion detection system. 

Accurate calibration can be achieved by means of a 
scheme in which an additional detector is used to count 
the electrons produced by ionisation |5fJ . If Nj atoms are 
ionised in a calibration, and the overall ion and electron 
detection efficiencies are rji and i] e respectively (where 
analogous to Eq. JTJ, r) e = r)det,ePr,e), then 



Ni = rnNi, 
N e = r] e Ni, 
N c = riiTjeNi, 



(5) 
(6) 
(7) 



where Ni is the number of ions counted, N e is the number 
of electrons counted and N c is the number of ion-electron 
coincidences counted. Therefore 



Vi 



Nc 
N K 



(8) 



which shows that only the coincidence and electron 
counts are required to find the ion detection efficiency; 
we don't need to know the efficiency of the electron de- 
tector. 

In order to count the number of coincidences, it is first 
necessary to specify how coincidences are identified. In 
general there will be a distribution of times taken for ions 
to move from the point of ionisation to striking the detec- 
tor and producing a pulse. The width of this distribution 
will depend on the range of initial velocities of the ions 
and the variation in electric field at different locations 
within the atom cloud. The same is true for the elec- 
trons. There will therefore be a distribution of intervals 



between ions and their corresponding electrons. With 
modern pulse counting electronics it is straightforward 
to produce lists of the arrival times of electrons and ions 
with precisions on the order of 100 ps [5l| and these can 
be used to generate a list of intervals between electrons 
and the following ions. If a histogram of these intervals is 
plotted, true coincidences will appear as a peak against 
a small background of false coincidences. This peak can 
then be used to define a window for true coincidences. 

The coincidence window may also contain false coinci- 
dences, which occur when an unpartnered ion (i.e. one 
for which the corresponding electron was not detected) 
arrives within the window of an unpartnered electron. 
Also, due to the finite size of the window, some true co- 
incidences may not be counted. The mean net number of 
false coincidences Nf and its uncertainty Of are derived 
in Appendix IbI 

In order to take false coincidences and background ions 
and electrons into account in the ion detection efficiency, 
we define corrected ion, electron and coincidence counts 
as NI = Ni ~ N ih , NI = N K — N Kb and NI = N r - 



Nf, respectively. Nn and N e b are respectively the mean 
numbers of background ions and electrons counted in the 
calibration period and can be measured by counting ions 
and electrons in the absence of the calibration atoms. 
The overall ion detection efficiency is then 



N' 

1 _ C 

n>: 



(9) 



The ion detection efficiency given by Eq. © is the value 
that applies in a calibration and is valid for the ionisation 
rate used in the calibration. To find the system ion de- 
tection efficiency (which is ionisation rate independent) 
we need to correct rji for the value of the ion temporal 
resolution probability p r that applies in the calibration. 
From Eqs and @ , the system ion detection efficiency 
is 



Vdet 



N' ePr 



and from Appendix O 



Pr 



l-T r Ni--{T r Nif, 



(10) 



(11) 



where r r is the minimum time interval for which two ions 
can be resolved from each other, and N is the average 
rate of ion detection in the calibration. 

We also need to determine the uncertainty in r/det- This 
uncertainty arises mainly from the fact that for a given 
number of atoms ionised, the number of ion, electron 
and coincidence counts are described by binomial dis- 
tributions. An additional contribution comes from the 
uncertainties in the number of background electrons and 
false coincidences counted during the calibration. Tak- 
ing all of these effects into account, and for values of p r 
close to unity, the uncertainty in the ion detection system 
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efficiency is 




(12) 



where a e b is the uncertainty in the number of background 
electrons counted and is equal to \J N e b, assuming the 
detection of background electrons is a Poisson process. 

To minimise the uncertainty in the ion detection ef- 
ficiency, it is desirable to minimise background counts. 
Background electrons and ions arise from three main po- 
tential sources. The vacuum gauges and ion getter pumps 
commonly used in ultra-cold atom experiments work by 
ionising atoms and hence some of the ionisation products 
may find their way to the detector. While the vacuum 
gauge can simply be turned off during the experiment, 
ion counts from the pump can be reduced by adding ap- 
propriate shielding. The third potential source of back- 
ground ions and electrons is the ionisation by stray laser 
light of the Rb atoms covering the interior of the vacuum 
chamber. 

Background ions and electrons can be prevented from 
reaching the detectors by using suitable ion optics, such 
as the apparatus presented in |48|. In this apparatus 
four electrodes are used to accelerate and focus ions and 
guide them to the detector. The ion trajectories were 
calculated using standard software and it was found that 
only a volume of ~ 1 cm 3 is imaged onto the detector. 
Ions produced outside of this volume are not detected 
and it is estimated that background counts can be almost 
entirely suppressed. 

From Appendix ]B\ the uncertainty in the number of 
false coincidences can also be made very small. From 
Eqs HI U|l and 1)12(1 . neglecting background counts and 
false coincidences, and for p r « 1, the average relative 
uncertainty in the ion detection efficiency is 



Vdet 



1 



Vdet 



VdetVdet,eNl 



(13) 



For calibration parameters Nj = 10 5 s _1 , rjdet = Vdet,e = 
0.8 and a calibration time of r ca ; = 10 s, this expression 
is accurate to within 1% for background electron and ion 
rates of up to 10 3 s _1 . Equation (|13fl shows that this 
calibration method works best for large atom numbers 
and high overall ion and electron detection efficiencies. 

In order to test the values for rjdet and a Vdct predicted 
by the calibration equations, numerical simulations of 
calibration experiments were carried out. These simu- 
lations included all the relevant features of a real calibra- 
tion, such as 

• ionisation following Poissionian statistics 

• detection of ions and electrons subject to their re- 
spective detection efficiencies 

• background electrons and ions 
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FIG. 3: Histogram of the inferred values of the ion detection 
system efficiency rjdet for 1000 runs of the calibration sim- 
ulation with rjdet = 0.9, electron detection system efficiency 
rjdet,e = 0.9, ionisation rate Ni — 10 4 s _1 , calibration time 
Teal = 10 s, pulse resolution time r r = 10 ns and, background 



electron and ion detection rates Nib 
tively. 



Net = 10 s respec- 



• resolution times of the detectors 

• generation of a coincidence window from the lists 
of electron and ion arrival times 

• use of the coincidence window to determine the 
number of coincidences 

Each simulation generated electron, ion and coinci- 
dence counts for both "real" and background electrons 
and ions. Then, as would be done in a real calibration 
experiment, the corrected counts were calculated using 
the estimated numbers of background electrons, back- 
ground ions and false coincidences. An inferred value 
and an inferred uncertainty of the ion detection system 
efficiency rjdet were then found using Eqs (fTTTf) and lJT2)l. 

Excellent agreement was found between the theory and 
simulation for a wide range of values. Figure |31 shows 
a histogram of the inferred ion detection efficiency for 
1000 runs of the simulation with rjdet = Vdet,e = 0.9, 
Ni = 10 4 s-\ T ca i = 10 s, r r = 10 ns and N ib = N, 



,b 



= W B" 



The inferred values of rjdet have a mean of 
0.899993 in excellent agreement with the actual value of 
0.9. The standard deviation of the inferred values was 
1.047 x 10^3, also in excellent agreement with the aver- 
age theoretical value of er I)dct , which for these parameters 
is 1.049 x 10~ 3 . 

To determine the values of rji and er, ;i that apply in 
a counting experiment, the calibrated values of rjdet and 
(Tri d et nrust be corrected for the value of the ion temporal 
resolution probability p r that applies for the ion detec- 
tion rate used in the experiment, r/i is found from Eq. 
and er r)i = cr I?det for p r as 1. For the case of ionisation 
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FIG. 4: Contour plot of the normalised count uncertainty n 
showing how k increases as the ionisation time tj decreases, 
for various combinations of initial values of the overall ion 
detection efficiency rji and its uncertainty a Vi . Crosses, circles 
and triangles correspond to ti — > oo, rj = 2 ms and ti = 1 
ms respectively. 



at a constant rate in atoms per second, p r is given by 
Eq. ijCTfl in Appendix IU1 If ionisation occurs at a con- 
stant per atom rate (i.e. ionising laser powers are held 
constant), we use the effective ion resolution probabil- 
ity p' r , given by Eq. I|C11|) . It should be noted that for 
the calibration to be valid, the ion collection efficiency 
rjcoii must remain constant between calibration and ex- 
periment. This could be achieved by calibrating using 
cold atoms in the same trap as used in the experiment. 
To achieve the lowest possible uncertainty in the atom 
count, rjcoii should be unity, which can in principle be 
achieved with well-designed ion optics. 



V. TRAP LOSS AND OPTIMAL IONISATION 
RATE 

In an atom counting experiment it is important to 
know how the rate of ionisation affects the normalised 
count uncertainty k. The greater the rate of ionisation, 
the smaller the value of p r and the larger its uncertainty. 
This in turn decreases rji and increases a Vi , resulting in 
an increased value of k. However in most experiments it 
would generally be desirable to ionise a cloud as quickly 
as possible without overly affecting the count accuracy. 
Figure 21 shows how k is affected by the ionisation time 
rj = Nj/Ni, for Nx — 10 3 and r r = 5 ns, assumed to be 
known with a 10% uncertainty. The crosses correspond to 
an infinitely long ionisation (in which case rji = rjdet) and 
the circles and triangles correspond to ionisation times of 
2 and 1 ms respectively. It can be seen that down to a 
certain ionisation time, k is close to its minimum value. 
However as the ionisation time decreases, k begins to in- 
crease rapidly. 



If there is negligible trap loss, the number of atoms 
ionised Ni is the same as the number of atoms originally 
in the trap N, and the inferred atom number Ni n f is the 
same as the inferred ion number Nj i n f. If there is trap 
loss however, we have to correct the inferred ion number 
to obtain an inferred atom number. In the remainder of 
this section we show how this is achieved and how the un- 
certainty in the inferred atom number is affected by trap 
loss. We also find the optimal ionisation rate which min- 
imises the uncertainty in the inferred atom number due 
to both trap loss and imperfect ion temporal resolution. 

To calculate the atom count corrected for trap loss, 
we assume a constant loss rate Ri and that ionisation 
occurs at a constant rate Ri (both rates are per atom). 
The respective time dependent probabilities of an atom 
being in the trap and of having been ionised follow 



P T = -(Ri+Ri)P T , 
Pi = RiPt- 



(14) 
(15) 



The probability that an atom is ionised rather than lost 
is 



Pi 



Pidt = 
i) Ri 



Ri 



Ri 



(16) 



The total number of atoms ionised Nj is described by a 
binomial distribution, which again can be approximated 
by a normal distribution with mean piN and standard 
deviation y /p I (l-p I )N H<|. The original number of 
atoms in the trap can then be inferred to be 



N, 



inf 



1+ j|] V '-/ 



with a statistical uncertainty of 



&N irlf ,stat 



(17) 



(18) 



Due to uncertainties in Ri, Rj and Ni t i n f, we get extra 
terms in the uncertainty in Ni n f and we find 
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' Ri 



N TMf , (19) 



with Ni t i n f and a^j inf given by Eqs © and ©. Thus 
in order to infer the atom number and its uncertainty, 
we need to know Ri and Ri and their uncertainties in 
addition to Ni, rji, and a rn . We now consider how these 
values can be measured in a given experiment. 

The sum of Ri and Ri can be found by plotting 
Ln[iVj — Ni(t)} against time, where Ni(t) is the cumula- 
tive ion count. The slope is then —(Rj + Ri). The value 
of Ri (and hence of Ri) may be found by plotting Rj + Ri 
against the intensity of the 778 nm laser for a number of 
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FIG. 5: The optimal 1/e time required to ionise the cloud 
is shown as a function of the atom number and the loss rate 
from the trap. 



runs. Ri is given by the extrapolated y-intercept of this 
plot. (Extrapolation is required because Ri cannot be 
measured by counting ions if no ionisation occurs. How- 
ever, Ri could be measured directly by absorption imag- 
ing at successive times in the absence of 778 nm light.) 

The optimal ionisation rate is that which gives the low- 
est value of k for given values of Ri, r)d e t, T r and N. In 
Fig. |5j the optimal value of 1/Ri is plotted as a function 
of TV for various values of 1/Ri, for rjdet = 0.9, <t., h = 0, 
T r = 10 ns and assuming a 10% uncertainty in the mea- 
sured values of Ri and Rj. The corresponding minimum 
values of n/n m in are plotted in Fig. [ft] where K. m in is the 
value of k for zero trap loss. It can be seen that reason- 
ably sized atom clouds (10 3 — 10 4 atoms) can be counted 
in times on the order of 10 ms with close to the theoretical 
minimum count uncertainties for realistic trap lifetimes. 

An improvement on the uncertainties resulting from 
ionising with constant per atom probability would be ob- 
tained by increasing the 778 nm laser power during the 
ionisation in such a way as to obtain a contant ionisation 
rate in atoms/s. This would mean that the cloud would 
be ionised faster with less trap loss. The time dependent 
rate equations describing this scheme could be integrated 
numerically to determine the probability that an atom is 
ionised rather than lost, from which the original atom 
number and its uncertainty could be inferred. 



VI. DISCUSSION 

We now consider the application of our atom counting 
scheme to the characterisation of unknown atomic sam- 
ples. To determine the number of atoms in a single cloud, 
one would simply count a certain number of ions Ni and 
use the known calibration data to obtain an inferred atom 
number and uncertainty Af m j ± 0jv <n , from Eqs © and 



FIG. 6: The minimum achievable value of the relative count 
uncertainty k expressed relative to K m in (the value for zero 
loss) is shown as a function of atom number and trap lifetime. 
For the given values of the overall ion detection efficiency r\i 
and its uncertainty a Vi , Kmin ~ 0.11. 

©. In the case that trap losses are not negligible, Eqs 
1|17JI and (|19|l are also required. If however we wish to 
know the atom number variance in a specific situation 
we need to analyse a number of repeated measurements. 
As an example, consider a sample of n measurements on 
atom clouds generated in an identical manner which, un- 
der ideal (noise free) circumstances, would yield identical 
numbers of atoms (for example usin g a q uantum tweezer 
to extract atoms from a BEC [32, [33| ). Each experi- 
ment yields a measure Ni n fj ± <JN inf ■ where j — 1, 2..n. 
The quantities of interest are the mean atom number, 
the atom number variance and the uncertainties in these 
values. 

The population mean can be estimated from the 
mean of all infe rred atom number s TV and has an 
uncertainty of s T /n + (a m N /-qi) 2 where sip is the 
bias-corrected sample total variance and is given by 
s T = l/(n — 1) Yl(Ninf,j — N) 2 . The true atom number 
variance a\ is estimated by the difference of the total 
variance of the measurements sp and the variance asso- 
ciated with the binomial statistics of the ion counting 
= N(l — rji)/i]i. Note that the expression for er^ 
does not include a term due to a r)i as in Eq. 0, since 
we assume that any error in the detection efficiency (cal- 
ibration) will be the same for each measurement. This 
leads to a systematic error in the atom count and does 
not contribute to the variance a 2 ^. 

The degree of number squeezing of the atomic popu- 
lation is described by the Fano factor Fa = a\/N p^ - 
A Poissonian sample has Fa = 1 and a sub-Poissonian 
(number-squeezed) sample has Fa < 1- We can define 
a total Fano factor Ft = s T /N and measurement Fano 
factor F m = er^/TV, which accounts for the fluctuations 
due to the binomial counting of ions. The population 
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Fano factor is then given by Fa — Ft — F mi and it can 
be shown that for a Vi < 5% and N > 5 the uncertainty 
in Fa is well approximated by op A = y/2/n{FA + F m ). 
Knowing F m from the calibration it is possible to de- 
termine how many measurements are required to obtain 
a suitably small uncertainty <jp A in the measured (nor- 
malised) atomic variance Fa- 

VII. CONCLUSION 

We have analysed photoionisation and ion detection 
as a means of achieving accurate counting of ultra-cold 
atoms. The scheme relies on a three-photon, cw ionisa- 
tion process and can also be used for efficient single atom 
detection with high spatio-temporal resolution. To count 
large numbers of atoms, one of the lasers acts as an opti- 
cal dipole trap for the target atoms. We have shown that 
a well calibrated ion detector with realistic efficiencies 
can achieve single shot atom counts with sub-Poissonian 
accuracies for atomic samples containing many thousands 
of atoms. This is potentially useful for studies of the 
quantum properties of cold atomic gases. The variance 
of an ensemble of systems can also be characterised in 
a manner that only weakly depends on uncertainties in 
the absolute detection efficiency. Similarly, this scheme 
could allow direct probing of atom number correlations 
in molecular down conversion and four wave mixing ex- 
periments which are predicted to show number difference 
squeezing. 

This work was supported by the Australian Research 
Council, Landesstiftung Baden- Wiirttemberg and the 
European Union (MRTN-CT-2003-505032). 

APPENDIX A: TABLE OF SYMBOLS 



N number of atoms originally in the trap 

N in f inferred no. of atoms originally in the trap 

ON inf uncertainty in the inferred atom number 

° N inf stat statistical uncertainty in the inferred atom no. 

Ni number of atoms ionised 

Ni.i n f inferred number of atoms ionised 

o-jv, inf uncert. in the inferred no. of atoms ionised 

Ni ionisation rate (s _1 ) 

r)i overall ion detection efficiency 

a Vi uncertainty in the overall ion detection eff. 

r) e overall electron detection system efficiency 

rj det ion detection system efficiency 

a t)det uncert. in the ion detection system eff. 

r/det,e electron detection system efficiency 

ri' dEt ion detector efficiency 

rjcoll ion collection efficiency 

■p r ion pulse temporal resolution probability 

p' r effective ion pulse temporal resolution prob. 

p r , e electron pulse temporal resolution probability 

k normalised count uncertainty 

Kmm normalised count uncertainty for zero trap loss 



T r 


pulse resolution time (s) 


Teal 


calibration time (s) 


Tl 


ionisation time (s) 


Twin 


coincidence window length (s) 


T 


interval between successive ion detections (s) 


Ni 


number of ions counted 


Ni{t) 


cumulative ion count 


Ni 


ion detection rate (s _1 ) 


n: 


ion detection rate as r r — > (s _1 ) 


Ni, C al 


ion detection rate in a calibration (s _1 ) 


N e 


number of electrons counted 


N c 


number of coincidences counted 


Nl 


corrected ion count 


K 


corrected electron count 


N' c 


corrected coincidence count 


N ib 


mean no. of background ions counted 


N eb 


mean no. of background electrons counted 




uncertainty in the background electron count 


Nf 


mean net no. of false coincidences counted 


°7 


uncertainty in the false coincidence count 


Pf 


probability of a given unpartnered ion being 




in the window of an unpartnered electron 


Pw 


probability of a true coincidence being outside 




the coincidence window 


<5int 


standard deviation of the time intervals be- 




tween electrons and their corresponding ions 


Ri 


trap loss rate per atom (s _1 ) 




uncertainty in Ri (s _1 ) 


Ri 


ionisation rate per atom (s _1 ) 


or.! 


uncertainty in Ri (s _1 ) 


Pt 


probability of an atom being in the trap as a 




function of time 


Pi 


probability of an atom being ionised as a func- 




tion of time 


Pi 


probability of an atom being ionised as t — > oo 


fint(T) 


probability density function of the time inter- 




val between successive detections 


/»(*) 


probability density function of arrival times of 




ions at the detector 





T r 7j det NRl 


n 


number of atom number measurements 


N 


mean of inferred atom numbers 


St 


bias-corrected sample atom number total 




variance 


Ni„f J 


inferred atom number in the j th measurement 




uncertainty in the inferred atom number in the 


j th measurement 


o\ 


atom number population variance 


_2 


variance in the atom counts due to the bino- 




mial statistics of the measurement 


Fa 


Fano factor of a population of atomic clouds 


<TF A 


uncertainty in the Fano factor of a population 




of atomic clouds 


Ft 


Fano factor of a sample of inferred atom 




numbers 


F m 


measurement Fano factor 
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APPENDIX B: NET NUMBER OF FALSE 
COINCIDENCES 

In this appendix we derive the required correction to 
a coincidence count to account for false coincidences and 
uncounted true coincidences. In a given calibration there 
will be Ni — N c unpartnered ions and N e — N c unpart- 
nered electrons. (We assume negligible background ion 
and electron numbers.) The probability of a given un- 
partnered ion being in the window of an unpartnered 
electron is well approximated by 

p f = (N e -N c )^, (Bl) 

T~cal 

where r win is the window length and r ca ; is the duration 
of the calibration. 

The probability of a true coincidence being outside the 
window is 

r-^te)' « B2 > 

where <ii n t is the the standard deviation of the time in- 
tervals between electrons and their corresponding ions. 

The net number of false coincidences counted in a given 
calibration is then 

N f = N f ±a f , (B3) 

where 

N f =p f (Ni-N c )-p w N c , (B4) 

and 

<Jf = yjp } (1 -pf) (Ni - N c ) +p w {l- Pw ) N c (B5) 

APPENDIX C: CALCULATION OF p r AND p' r 

As the overall ion detection efficiency rji is ionisation 
rate dependent, relating the ion detection efficiency mea- 
sured in a calibration to the ion detection efficiency which 
applies in an experiment requires finding an expression 
for p r in terms of the ionisation rate [see Eq. JJ)]. We 
consider two cases: constant ionisation rate in atoms per 
second, and constant per atom ionisation rate. We look 
first at the former. 

We assume the ion detection system has some resolu- 
tion time r r , such that if an ion arrives at the detector 
within r r of the previous ion it will not be counted. We 
define N* as the average rate (ions/s) at which ions are 
detected as r r — > 0. Assuming a constant average de- 
tection rate, the detections follow Poissonian statistics, 
and the time interval between successive detections r, is 
described by a negative exponential probability density 
function [49| : 

/int(r) = N*e^' T for r > 0. (CI) 



The probability that an ion arriving at the detector is 
resolvable from the previous ion is 

p r = l- f T fint(r)dT. (C2) 
Jo 

From Eqs (□} and (JS2) we find 

p r = exp(-T r N*). (C3) 

The rate at which ions are actually detected is 

Ni=p r N*. (C4) 

The number of ions counted during a calibration period 
T C ai is described by a Poisson distribution: Ni = NiT ca i ± 
\J NiT ca i ■ If Ni ions are counted in time r ca ;, the ion 
detection rate in the calibration is inferred to be 

N itCal = ± ^5. (C5) 

From Eqs JU3J) and (EH) we find 

N* = (1 + T r Ni)Ni + O ((r,.^) 3 ) ■ (C6) 

Substituting Eq. I|C6|I into Eq. (|C3H and neglecting the 
third order term (which results in a relative error in p r 
of less than 10 -9 for r r Ni < 10~ 3 ) we get 

p r = l-T r N i -^(T r N i ) 2 . (C7) 

Now we consider the case of ionisation with a constant 
per atom ionisation rate. In this case we need to use an 
effective ion resolution probability p' r because by ionising 
with a constant per atom rate, the ionisation rate in ions 
per second is a function of time and hence so is p r . The 
effective ion resolution probability is given by 

/>oo 

Pr = / Pr(t)fa(t)dt, (C8) 
JO 

where from Eqs [|C3f) . ((TIJ) and (JTSJ 

p r (t) = exp (-/?e-^+^*) , (C9) 

where (3 = T r r)detN Ri and 

f a (t) = (Ri + R l )e^ R ' +R ^ t (CIO) 

is the probability density function of arrival times of ions 
at the detector. It can be shown that, neglecting terms 
ofO(/3 3 ), 

Pr = l-\P+\P 2 - (Cll) 
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